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Abstract 

We perform an OPE analysis of the flavorless non-perturbative gluon propagator and 
the symmetric three-gluon vertex in the Landau gauge. The first subdominant operator is 
A^A^ which can condensate in the Landau gauge "vacuum" although being a non-gauge 
invariant operator. We neglect all higher dimension operators. Then the gluon propagator 
and the symmetric three gluon vertex only depend on one common unknown condensate. 
We propose a consistency check from lattice data. At two loops for the leading coefficient 
and with 1/p 2 corrections at tree-level order the two fitted values for the condensate do 
not agree. At three loops we argue that the today unknown (3% lOM should be equal to 
f .5(3) x l3f OM = 7400(f500) to fulfill the OPE relation. Inclusion of the power corrections' 
anomalous dimensions should improve further the agreement. We show that these techniques 
cannot be applied to the asymmetric three gluon vertex with one vanishing momentum. 

P.A.C.S.: I2.38.Aw; I2.38.Gc; I2.38.Cy; 11.15.H 

1 Introduction 



Much effort has being devoted in the last years to the non-perturbative computation of the 
running QCD coupling constant a s (p) A most surprising result arising from the lattice 

computation of the gluon three-point Green function [[jj is that a ~ 1/p 2 power correction is 
still sizeable for momenta of the order of ~ 10 GeV !! 

These power corrections were qualitatively discussed in |Q] in connection with several physics 
phenomena such as interquark potential, the gluon condensate puzzle, Landau pole, etc... It was 
also strongly stressed that a s {p) being non-perturbatively denned on a lattice in the particular 
Landau gauge, a contribution to the operator product expansion (OPE) from a gauge dependent 
local operator should not a priori be neglected. More precisely, the A 2 = A^A 11 operator can 
acquire a vacuum expectation value (v.e.v.) because of the non-gauge-invariant nature of the 
vacuum. As a matter of fact, the authors of ref. || showed in a general way that the OPE 
for QCD propagators unequivocally generates non null v.e.v. of A^A^. The important point is 
that in contrast to that of the gauge invariant G^G^y generating 1/p 4 powers, it produces 1/p 2 
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corrections which may be of comparable size to higher perturbative correction (three-loops, for 
instance) for the considered momenta. 

The key remark here is that A 2 is the only dimension-two operator able to generate a non 
vanishing v.e.v. The purpose of the present paper is to exploit this fact according to OPE 
philosophy in the following manner: we will study both the gluon propagator and the three 
gluon vertex (two- and three-point Green functions) in order to examine whether the same (A 2 ) 
condensate can consistently describe the evaluated power corrections to both quantities. For 
rather obvious reasons which will be explicited later, we will totally neglect 0(1/ p 4 ) corrections. 

The OPE approach consists in expanding the Green functions in inverse power series in p 2 , 
every term corresponding to a local operator. The coefficients which multiply the operators are 
calculable in perturbation In practice of course the coefficients are expanded to a given finite 
order, not only because the perturbative series are only asymptotic, but also because computing 
higher orders is a rapidly increasing task which has not yet been performed. 

The potential dangers when using OPE to relate different Green functions via the same 
condensates has been stressed by several authors [11, 1^]: it was claimed that the coefficients of 
the leading operator at any given momentum scale have to be computed to a sufficiently high 
perturbative order to make the perturbative uncertainty smaller than the power corrections. 

Our input data are lattice computed Green functions [^-0], covering a momentum range from 
2 to 10 GeV. We will analyze them according to OPE. The gluon propagator and the symmetric 
three-point Green function (p 2 = p\ = p 2 ) can be expanded into perturbatively computable 
coefficients and the only unknown (A 2 ) condensate. This allows to perform a test of OPE in 
a relatively favorable situation. Our first "ace" is indeed this wide momentum window which 
reduces the risk of confusion between neglected perturbative (logarithmic) contributions and the 
~ 1/p 2 termsH. Our second "ace" is that 1/p 2 terms are rather easy to display in the 2-10 
GeV window |7j while, for example, 1/p 4 contributions would only influence the lower part of 
this window making the fits hazardous. Our third "ace" is the high accuracy achieved by our 
lattice calculations. 

The leading coefficient (the coefficient of identity operator in the OPE) of the gluon prop- 
agator is known up to three loops but the symmetric three-gluon vertex is only known at two 
loops. To ask for an OPE consistent description when leading power coefficients are only ex- 
panded at two loops might be exposed to the danger mentioned above [12|. We intend to check 
this. We will also try the three-loop consistent description by fitting the unknown perturbative 
information for the three-gluon vertex, the @2 parameter of beta function in MOM scheme. The 
comparison with future perturbative computations will obviously be in order. 

We have only computed at tree-level order the perturbative coefficients multiplying the A 2 
operator. We will therefore perform the fits with a constant times 1/p 2 term. A knowledge of 
the anomalous dimension of the coefficient which multiplies 1/p 2 will be welcome [Jllfl. 

The theoretical setting of our use of OPE is described in Section ^. A preliminary analysis 
of results from previous lattice data are presented in Section ||. We finally discuss and conclude 
in Section 



2 Still it should be noted that the recent perturbative computation of /?3 in asymmetric MOM scheme[f4|] leads 
to a four-loop estimate of the empiric power coefficient introduced in ref. jjj wich is roughly 0.8 times the three- 
loop one. This shows that some amount of disguise of higher order perturbative terms into power corrections is 
difficult to fully avoid. 
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2 OPE for the gluon propagator and a s (p) 



The purpose of the present section is to develop an OPE-based model describing the non- 
perturbative power behaviour of gluon propagator and a s (p). To this goal, we will compute 
the first power correction in OPE (always at tree-level) for two and three gluon operators, by 
following the prescriptions of [p7Cf| . Some general results for the gluon propagator in full QCD can 
be found in ||. For completeness we will describe our own calculations for both the propagator 
and the three gluon vertex. We have checked agreement for the propagator with M. 

The two- and three-point bare lattice Green functions are renormalized by imposing mo- 
mentum subtraction prescriptions. Let us insist: such a renormalization prescriptions leads 
to the definition of a gauge-dependent coupling which consequently admits a non vanishing 
contribution from the non-gauge invariant OPE A^A^. 

We work in the pure Yang-Mills QCD, without quarks. OPE yields 



T <(-p)4(p) 



(cb)£(p) 1 + (ci)2£(p) : <'(0) : + (c 2 ) a u Sb) : Aj(0) 4',(0) 



fj.ua' b' 



+ 



(1) 



r(l2(pi)^(p2K(P3)) = 

(do)^ P (Pi, M ) 1 + (di)^(pi,p a ,ft,) : <,(0) : 

+ (d2fXi(Pi,P2,P3) :A£(0)j#(0): +... ; 

(2) 

written only in terms of normal products of local gluon field operators, where A (A) stands for 
the gluon field in configuration (momentum) space, a, b being colour indices and fj,, v Lorentz 
ones. The notation T() simply refers to the standard T* product in momentum space. The 
normal product of Eqs. will be defined below in reference to the perturbative vacuum. 

The r.h.s. in eqs. (0,0) contains operators of dimension up to two. We have omitted local 
operators including, for instance, gluon field derivatives or ghosts. Indeed at the order 1/p 2 , 
which we consider in this letter, they have a vanishing v.e.v. (Obviously < d^A^ 1 >= 0). 



2.1 Gluon propagator 

The following step is to take the QCD v.e.v. (in our flavorless universe) in both r.h.s. and l.h.s. 
of equation ([!]), 



(t (A«(-p)A b u ( P ))) = (co)f u (P 2 ) + (c 2 ) a Xr b : * a ' 6 Vv 4( ^f 1 !) +••• ; (3) 

where A 2 = : A";(0) A^, (0) : . A first immediate simplification in r.h.s. of Eq. ([|) with 
regard to Eq. (|l]) is the disappearance of terms containing an odd number of local gluon field 
operators. They contribute only at higher orders in the OPE because of the impossibility to 
build a Lorentz invariant tensor (with odd number of indices) without any external non-null 
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momentum. The normal product in eqs. ([[],§) refers to the perturbative vacuum |0 >. To be 
more specific the non-perturbative condensate < A 2 > is defined from eq. (^) as follows ^ : 
given the l.h.s computed non-perturbatively on a lattice?], given the coefficients cq and C2 in 
the r.h.s computed perturbatively at a given order in perturbation theory and for a particular 
renormalization momentum scale and scheme < A 2 > is unambiguously defined. 

The standard method to compute the perturbative expansion of the OPE Wilson coef- 
ficient is to compute to the wanted order an appropriate matrix element of the l.h.s. operator 
of Eqs. (jl],§). The perturbative gluon propagator 

(G { S rt )t(v 2 ) = (0\T (Al(-p)A b u {p)) |0> (4) 
at the considered order in perturbation theory yields c$ while C2 may be obtained from 

(g\T (A«(-p)A b u (p)) ^connected = (c 2 )J£(p) (g \ I A}{0) j#(0) : \g) (5) 
with two soft external gluons. From 

(g\A 2 \g) = 2 + 0(a), (6) 

we may compute 02- At tree-level order the 0(a) is neglected and eq. (|5|) allowing to compute 
C2 from the computation of the l.h.s. : ^] 



, M 1 {A^)Al{-p)A b v {p)A%m 

The ratio in the r.h.s of Eq. (|7|) represents symbolically diagrams with four gluon legs where the 
two of them carrying zero momentum are explicitly cut. Consistently with Eq. (|6|), we compute 
in Landau gauge the tree-level diagrams shown in fig. [l| to obtain and then, for p 2 = —k 2 , 

we can write: 



(T (a% 



-P)%<P))) = -zh^- P -^)5^ 



.2^(2) fu 2, , 3g 2 (A 2 ) ! , n / _ +1 „ A 2 



As already mentioned, there is a series of problems in the practical use of OPE tried in 



this paper [11], [12]. In general, the whole expansion in terms of condensates suffers from the so- 
called renormalon ambiguities : the resummation of the non-Borel-summable perturbative series 
for the leading operator Wilson coefficient is ambiguous, this ambiguity being compensated by 
the ultraviolet renormalon ambiguity of the matrix element of the subdominant operator A 2 . 
By luck the infrared renormalon ambiguity cancels off when comparing the contribution of the 
condensate < A 2 > in both the propagator and the three gluon vertex. In this way we may 

3 We thank Y. Dokshitzer and G. Korchemsky for illuminating discussions on this point. 
4 Being a gauge dependent quantity it cannot be accessed from experiment. 

5 Beyond tree-level order the radiative corrections, 0(a) terms in (H), are calculable, but we will skip this issue 
in this paper. 
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Figure 1: Four gluons tree-level diagrams contributing (with all their possible permutations) to the 
(c2) a ^J t ° OPE perturbative coefficient of the gluon propagator. Crosses mark the gluon legs due to the 
external soft gluons. 

simply decide to truncate the perturbative series to some order and the universality of the 
condensates is preserved provided this order is clearly stated. 

A remnant of this renormalon problem arises from the truncation of OPE and perturbative 
series for the Wilson coefficients. One has to be sure to have gone far enough in perturbative 
expansions and in the inverse power expansion to have a coherent scheme. We will try to justify 
later our truncation of the power series to order 1/k 2 . For the perturbative series of the leading 
contribution we have to choose a truncation order n according to a compromise: it should be 
high enough to make sure of the dominance of non-perturbative contributions compared to the 
perturbative uncertainty from neglected higher orders; but low enough not to destroy the good 
asymptotic behaviour of the series. 

The gluon propagator turns out to be known up to the third loop in its perturbative expansion 
PkpHfl, but not yet the symmetric MOM coupling constant. We can thus try the exercise with 
both n = 2 and n = 3, but making in the latter case some fit for the unknown three loop 
coefficient (/?2)- This is not too high an orderj^J over the considered momentum window but 
it might be too low. Concerning power corrections higher than 1/k 2 : our fits do not show the 
necessity for 1 /k terms (which would exhibit a very distinguishable functional dependence) over 
the large momenta window we consider. This is lucky since in practice it would be hopeless to 
deal with 1/k 4 terms, too many operators being needed as well as, for consistency, many higher 
orders in the perturbative expansion of the Co, do, ci and di coefficients. Anyhow the consistency 
of the matching to lattice non-perturbative data for the two- and three-point Green functions 
will allow to test the validity of our assumptions. 

In brief, within the scope of the present paper we will satisfy ourselves by modeling power 
contributions to our Euclidean two-point Green function, 

(A%{-p) At(p)) = 5 ab (V - ^) G^(p 2 ) , (9) 

such that 

6 In the MOM scheme the fourth loop coefficient of the /3 function is known[[L4j and the series looks like being 
in its perturbative regime for momenta higher than ~ 3 GeV 
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(10) 



where Z^^^ = k 2 G n i oops , n = 2, 3 and by studying systematically the stability of the matching 
of lattice results to formula ( JlOj) for different choices of the momentum window. 

All the quantities in Eqs (|3}JlO|) are bare objects, where the explicit dependence on the 
regularization cut-off has been omitted for the sake of simplicity. However, applying MOM 
prescription at a certain renormalization scale // 2 , a renormalized version of Eq. ( |10| ) is obtained, 



tfG^fcV) = « PS (^ 2 ) + c ^ 2 ) ^ri) > (u) 

where, in MOM philosophy, (A 2 ) R/1 2 is renormalized at fi 2 such that the Wilson coefficient 
takes the tree-level value computed from the diagrams of fig. [j] : c {l^i^) = 3(<?/?(/i 2 )) 2 / '/x 2 . 
■^R f nioops(^ 2 ' m ^1- ( PI ) behaves according to the three-loop (two- loop) anomalous dimension 

of the gluon propagator Q and is such that ^ 2 G^(fi 2 , /i 2 ) = 1. This last condition imposes 
an overall normalization which is totally irrelevant for our lattice data fitting strategy. In this 
respect the situation will slightly differ in the next section. By definition the running k 2 is 
there the renormalization point for the running coupling constant. It is then the condensate 
renormalized at k 2 , {A 2 ) Rk 2, and the coefficient accordingly renormalised at k 2 which have to 
be dealt with. This has to be taken into account when computing the anomalous dimension 
of the factor multiplying 1/k 2 in each case as will be done soon [p.5[j. Nevertheless, the goal 
of the present paper is to develop the principle of an OPE testing procedure and to perform a 
preliminary check. In the scope of this preliminary check, where only tree-level power corrections 
are taken into account, we will simply take zero anomalous dimensions : 

c(k 2 ^ 2 )(A 2 ) R ^=3g 2 R { -^, (12) 
where g R (A 2 ) R is taken to be a non-perturbative constant which can be fitted from lattice data. 
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2.2 Running coupling constant 

In momentum subtraction schemes, in Landau gauge, the renormalized coupling constant can 
be defined at the renormalization momentum scale k 2 as (See for instance ||) : 



G^\p 2 ,plp 2 )(Z M OM {k 2 ))3 ^ ^ 



G^{p 2 )G^{p 2 )G( 2 ){p. 



where Z MOM and G^ were defined in the previous subsection, and G^ (p 2 ,P2,p 2 ) is the three- 
point scalar form factor. Eq. (|l^) is directly written for an Euclidean metric. Nevertheless, 
as done in the previous section, we prefer to compute the scalar form factor in the more famil- 
iar Minkowskian metrics and translate the final result to Euclidean. This form factor should 
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be obtained by projecting of the three-point Green function, {A a ^{pi)A h l ,(p2)A c p (p^)), over the 
transverse tree-level three-gluon tensor, 



/rptree\ 
V 1 



(pj _ p2 y' 3 + cycl.perm. [J U 



i=l,3 



'Mi ft 



(14) 



considered for a given configuration of the three momenta with an overall renormalization scale 
k 2 . In the following we shall consider the two standard momentum subtraction kinematics: 
= p\ = pi = -k 2 (MOM) and p\ = p\ = -k 2 , p\ = (MOM). 
In the MOM case, the scalar form factor, (k 2 , k 2 , k 2 ) , is obtained from: 



Pi 



G^(k 2 ,k 2 ,k 2 ) f 



abc 



1 



(A a J Pl )Al(p 2 )A c Jp 3 )) 



tree\^i/i2/i3 



/ rptvee^ 



+ 



-18k 2 " " 

(Pi ~ P2) P (P2 ~ P?,Y{P?, ~ PlY 

-2k 2 



(15) 



Then, we follow from Eq. (g) a procedure to obtain the OPE perturbative coefficient cfo anal- 
ogous to that to the one previously used for the gluon propagator. The diagrams to compute 
are those plotted in Fig. ||. In total analogy to the gluon propagator case, the result for the 
tree-level OPE correction of the three-point Green function, translated to Euclidean metrics, 
leads to 



G^(k 2 , k 2 , 



k 2 ) 



G 



(3) 
pert 



+ 



4(iV2-l) k 2 



(16) 



where, 



again, GpJ Tt stands for the purely perturbative three-point form factor. 
Then, we insert Eqs. (10[l^) in the MOM coupling constant definition, Eq. (|l3|), where the 
Green functions scalar form factors are automatically renormalized in MOM prescription at the 
scale k 2 by the 3/2-power of z M , and finally obtainf] 



(3) 



a. 



MOM 



(k 2 ) = (a s MOM (k 2 )) 



pert 



1 + 



(17) 



One interesting remark is that all the OPE contribution to a s MOM (k 2 ) defined by Eq. ( |i~3| ) is 
coming from the explicit renormalization factor (2 ,MOM )3/2_ After amputation of the external 
legs, diagrams (a) and (b) in fig ^ are obviously canceled, and a direct cancellation of the 
irreducible diagrams (c) and (d) happens. 

As already mentioned only^j the two-loop coefficient of the beta function is known in the 
symmetric MOM scheme. The gluon propagator anomalous dimension is known at three loops 



7 The 0(a)-corrections in the second r.h.s. of Eq. ([17 
in our tree- level OPE correction. 

8 the two first coefficients, /3q and /3i are universal 



are of the order of the assumed theoretical uncertainty 
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Figure 2: Five gluons diagrams to be computed, among all the possible new ones obtained by permutations 
of leg labels, in order to determine the tree-level OPE correction to the gluon three-point Green function. 



and previous works || has revealed the important role played by the third loop in its appro- 
priate perturbative description. Moreover, a previous analysis J5J including power corrections 

of a s MOM (k), in MOM scheme, also stressed the necessity of working at least at three-loops to 
get a good estimate of This is why we shall add to the two loop test of OPE a three-loop 

one with a fitted fo- The latter test is performed consistently at the three- loop level both for 
the gluon propagator renormalization constant, Z MOM (k 2 ), and for the perturbative part of 
a s MOM (k 2 ), which we will take to be the numerical inversion of 



A- 



k 



MS 3.334 



exp 



-2vr 



x exp 



P oa MOM 

10k Va 



( Poa 
1 4vr 



MOM 



arctan 



2ft +/3 2 a s MOM /47r 



arctan 



MOM 2 " 



A 




^5 



(18) 



where A = 2/?o/?2 — 4/3f , the unknown coefficient being fitted. The interplay between the 
ignorance of non-tree-level OPE correction and this fitted three-loop information is hoped to be 
controlled by studying the stability of the fits for different momentum windows. 

The OPE approach used for propagator and symmetric three-point Green function is based 
on the factorization of the non-perturbative soft gluon contribution (related to the vacuum 
properties) and the perturbative hard gluons in the external legs of the Green functions. In 

the MOM case, things are a little different. That has to be adapted since one of the external 
gluons of the Green function carries zero momentum. OPE large momentum expansion can 
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nevertheless be done for the two other gluon fields carrying non-null momentum by following 
Eq. ([!]). As only an odd number of gluon field operators permits to build a Lorentz invariant 
v.e.v., non-null contributions will come now for terms of the Wilson expansion with an even 
number of local operators. Then, we will have 

(T (A*{-p) A b v {p)A^))) = (c x )%{p 2 ) (: <,(0) : ^(0)) + . . . . (19) 

As for the propagator or the symmetric vertex, it is immediate to prove that (ci)**^(p 2 ) takes the 
tree-level expression for the asymmetric three-point Green function (with the explicit amputation 
of the zero- momentum gluon leg). Concerning subleading OPE corrections, they should be 
obtained from local operators of higher dimensions, and then a "zoo" of different condensates 
appear. In principle, ( : A^A^A^, :A c p ),(: d^c a 'c b ' : A°) or ( : d p ,A J (0)<, (0) : A°), are all 
candidates to give non-null v.e.v. contributions. Furthermore, the problem due to the soft gluon 
in the vertex definition re-appears now in a new way: no v.e.v in the r.h.s of Eq. ( |1S| ) involves 
directly the gluon condensate (A 2 ). Of course, this does not forbid an explanation of power 

corrections in terms of OPE for asymmetric a s MOM , but their non-perturbative condensates 
cannot be directly put in relation with those from symmetric a s MOM and gluon propagator. 
Consequently, with no additional phenomenological assumption no consistent test shall come 
from the running coupling constant in the asymmetric MOM scheme and it will be left aside in 
our present analysis. 

3 Fitting data to our ansatze 

Our goal is now to try a consistent description of lattice data for two- and three-point Green 
functions introduced in ref. |5]-|7| from the previously developed formulas for power corrections. 
Details of the lattice simulations, of the elaborated procedures to obtain an artifact-safe data 
set or of the definition of regularization-independent objects permitting lattice regularized data 
to be matched to continuum quantities in any scheme, can be found in those references. We 
directly show here the results of the matchings. 

In a previous work jjj, a fit of both the empiric coefficient for a power contribution and A^jg 

were simultaneously performed by matching lattice data for a s MOM (k 2 ) to a formula including 
non-perturbative additive power corrections. The A^jg parameter thus obtained, A^jg = 2371 13 
MeV, is in total agreement with the one estimated by ALPHA group, A^jg = 238(19) MeV. The 

same kind of fit of a s MOM (k 2 ) can be done by assuming the power correction to be proportional to 

a s M (fc 2 ) analogously to what is done in the MOM case in eq. ([T7| ). The fitted Aj^g parameter 

raises in the latter fit to ~ 250 MeV. This difference coming from the assumed oc a s M M (k 2 ) 
dependence of the 1/k 2 term tells us something about a certain additional systematic uncertainty, 
not taken into account in ref. (?|, due to neglecting anomalous dimensions of power corrections. 

We will take the central value of the estimated A-^g in ref. 0, A^g = 237^^ MeV, for 
a MOM an( j gj uon propagator. 



Two loop fits We have fitted the gluon propagator and the symmetric coupling constant 
according to eqs. ( pT|p^ ) and ( |i"7| ) with the leading perturbative parts computed to two loops. 
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In both cases we fit the condensate. The result is: 



propagator: g 2 R (A 2 ) R = (2.32(6)GeV) 2 ; a MOM : g 2 R (A 2 ) R = (4.36(12)GeV) 2 . (20) 



with x 2 /d.o.f ~ 1.2. The clear disagreement between both estimates could be a sign that we 



are in the situation indicated in [12] that the perturbative order is too low to give an acceptable 
precision in the estimate of the power corrections. We need to check if going to three loops is 
enough. 



pMOM is not k nown _ Then, our strategy will be to extract the value for the 



Three loops fit ^ 2 

gluon condensate, (A 2 ), from the matching of ( pl|Jl2]) to gluon propagator lattice data; and then, 
by using that value, to search whether any sensible value of the unknown MOM coefficient 02 
allows to reproduce data for a s M (p). 

In fig. | we plot the best fit of eq. (0®> for 



g 2 R (A 2 ) R = (2.76(4)GeV) 2 



(21) 



to gluon propagator data. This error is computed by jacknife's method. The x 2 is 0.86 per d.o.f. 
in that fit, slightly better than for the two- loop fit. It is amusing that the three-loop condensate 
stemming from the propagator is not so different from the two-loop one. 
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Figure 3: Formula with the assumed value of (237 MeV), fits impressively the data for g 2 (A 2 ) = 
(2.76(4)GeV r ) 2 over a large momenta window from 2.6 GeV to ~ 10 GeV (fig. a). The optimal window is 
fixed by requiring a reasonable stability of the fitted parameter (g 2 (A 2 ) ) for different choices of the lower 
momentum (fig. b). 

With this value of the gluon condensate we find that 

pMOM = L53(?) x ^MOM _ 740Q ± 30 Q (22) 

optimises the quality matching (x 2 /d.o.f. = 1.1) of formula ( |i~7| ) to lattice data (see fig. ^). 
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Figure 4: A^jg- assumed to be 237 MeV and g 2 (A 2 ) — 0.236(7) from gluon propagator analysis, the best 

fit of Eq. (H|) to lattice data (gives f)f OM = 1.53 x f3f OM (fig. a). Fig. b shows the stability of this 
parameter faced to different choices of the lower momentum of the fitting window. 



4 Discussion and conclusions 

A certain underestimation of the errors presented in the previous section should be admitted as 
a consequence of fixing Aj^g to its central value 237 MeV of [[7j without errors. We can correct 
for this by the same fitting strategy applied with different starting assumptions for Aj^g, within 
its error interval, this results in a a wider dispersion: 

pMOM i pMOM = L5(3)> (23) 

Ajyjg and fi^ iOM appear to be strongly correlated. The goal of the present paper is less to give a 
reliable prediction of a certain perturbative parameter, as /?^ OM , than to carry out a preliminary 
study of the inner consistency from a description in terms of OPE predictions of lattice data. 
Let us repeat, the knowledge of three- loop information in the MOM scheme and the anomalous 
dimensions of the power corrections are prerequisite to any firm conclusion^. The comparison of 
our estimate of /3^ OM and the future perturbative computation of this parameter^] are totally 
necessary to conclude whether signs of such a consistency are present or not. 

Our analysis of lattice data from gluon propagator and running coupling constant unequiv- 
ocally proves power corrections (1/p 2 ) to be necessary to match the data. The dimension 2 of 
power corrections clearly manifests the presence of non-gauge-invariant condensates (in partic- 
ular, in our pure gauge QCD analysis, A^A^), if OPE is invoked to explain the results. Since 
we work in the Landau gauge the appearance of such a non-gauge-invariant condensate is in full 
agreement with all theoretical expectations. The value for ygR < A 2 >r around 3 GeV might 
seem unexpectedly large but this is a subjective matter. The perturbative Wilson coefficients 

9 A very recent lattice calculation within the two-dimensional non-linear sigma model |TJ| insists on the necessity 
of knowing the OPE coefficients at least to one loop to reach agreement with OPE predictions. 
10 A computation of may appear soon J17[. 
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for the OPE of the involved Green functions have been theoretically computed. Then, in order 
to satisfy the universality of the non-perturbative gluon condensate appearing directly in OPE 
of both gluon propagator and MOM coupling constant, we need to assume /?^ OM = 7400(1500). 
The perturbative computation of this parameter as well as that of the anomalous dimensions of 
the power corrections will either exclude the interpretation of power corrections just in terms of 
universal gluon condensates at three loops, or on the contrary throw a positive sign supporting 
the usual U QCD sum rules' 1 procedure to connect non-perturbative effects for different correla- 
tion functions. All the theoretical tools for such a OPE testing procedure have been developed 
in the present paper. 
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